Neutrino spectra evolution during proto-neutron star deleptonization 
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The neutrino-driven wind, which occurs after the onset of a core-collapse supernova explosion, 
has long been considered as the possible site for the synthesis of heavy r-process elements in the 
Universe. Only recently, it has been possible to simulate supernova explosions up to ~ 10 seconds, 
based on three-flavor Boltzmann neutrino transport. These simulations show that the neutrino 
luminosities and spectra of all flavors are very similar and their difference even decreases during 
the deleptonization of the proto-neutron star. As a consequence, the ejecta are always proton rich 
which rules out the possible production of heavy r-process elements (Z > 56). We perform a de- 
tailed analysis of the different weak processes that determine the neutrino spectra. Non-electron 
flavor (anti)neutrinos are produced and interact only via neutral-current processes, while electron 
(anti)neutrinos have additional contributions from charge-current processes. The latter are domi- 
nated by ^-absorption on neutrons and i/ e -absorption on protons. At early times, charge-current 
processes are responsible for spectral differences between v e , z? e and ^ M / T . However, as the region 
of neutrino decoupling moves to higher densities during deleptonization, charge-current reactions 
are suppressed by final state Pauli-blocking. u e absorption on protons is suppressed due to the 
continuously increasing chemical potential of the neutrons, absorption on neutrons is blocked 
by the increasing degeneracy of the electrons. These effects result in negligible contributions from 
charge-current reactions on timescales on the order of tens of seconds, depending on the progenitor 
star. Hence, the neutrino spectra are mainly determined from neutral-current processes which do 
not distinguish between the different flavors and results in the convergence of the spectra. These 
findings are independent of the charge-current reaction rates used. It rules out the possibility of 
neutron-rich ejecta at late times and the production of heavy r-process elements from non-rotating 
and not magnetized proto-neutron stars. 

PACS numbers: 26.30.Jk 97.60.Bw 



I. INTRODUCTION 

Core collapse supernova explosions of stars more mas- 
sive than 8 M Q is an active subject of research in theo- 
retical astrophysics [l[ . They are related to the revival of 
the stalled shock, which forms when the collapsing core 
bounces back at nuclear densities. Right after the forma- 
tion, the shock looses energy due to the dissociation of 
heavy nuclei and to the release of the i/ e -deleptonization 
burst produced by electron captures on protons. At this 
moment the shock becomes an accretion front and stalls 
at a radius between 100 and 200 km. Different mech- 
anism have been suggested to revive the shock and fi- 
nally explode the star. These are the magneto-rotational 
mechanism the acoustic mechanism due to the 

dumping of sound waves § and the delayed neutrino- 
heating mechanism [f| . Recently, it has been shown that 
a phase transition to deconfined quark matter can lead 
to an additional collapse during which a second hydrody- 
namic shock wave forms that triggers the explosion @,@]- 

Independently of the explosion mechanism, a proto- 
neutron star (PNS) forms at the supernova core immedi- 
ately after core bounce. The initially hot and lepton-rich 
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PNS cools via continuous emission of neutrinos of all fla- 
vors in a period of several tens of seconds after the explo- 
sion has been launched. Neutrino-absorption processes 
at the PNS surface deposit energy, which drives a matter 
outflow known as the neutrino-driven wind. It has been 
investigated as a possible site for the production of heavy 
r-process elements | g.[l0| . These pioneering works where 
followed by analytic parametric [l3j and steady state 
wind models [lj, [l5[ . They showed that the production 
of heavy r-process elements requires a mass outflow on 
a short dynamical time scale (a few milliseconds), high 
entropies per baryon (above 150 ke) and a low proton- 
to-baryon ratio, i.e. electron fraction Y e < 0.5. Recent 
hydrodynamics simulations [H|,[l3] showed that short dy- 
namical timescales can in fact be achieved but fail how- 
ever to obtain the necessary entropies at times relevant 
for r-process nucleosynthesis [l8j . 

The nucleosynthesis outcome is very sensitive to the 
electron fraction of the ejected matter [2(| that is de- 
termined by the competition between electron neutrino 
absorption on neutrons and antineutrino absorption on 
protons and their inverse reactions. Deep in the PNS in- 
terior, neutrinos are in chemical equilibrium with matter. 
However, at the PNS surface where temperature and den- 
sity drop, neutrinos decouple from matter. Since \i and 
r neutrinos interact only via neutral-current reactions, 
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they are the first to decouple at the highest tempera- 
ture. For the very neutron-rich conditions found at the 
PNS surface, 9 e decouple before v e . The neutrino spec- 
tra reflect local properties of matter at the position where 
they decouple. Hence, one expects the follo wing hierar- 
chy of neutrino energies: £ v > £p c > £ ^e [2-l|> where 
e = (E 2 )/(E). (E) is the mean neutrino energy and 
{E 2 ) is the square value of the root- mean-square (rms) 
energy. The fact that v e have larger average energies 
than i/ e suggest that the ejecta will be neutron rich. How- 
ever, one also has to consider the other energy scale in 
the problem, i.e. the neutron-to-proton mass difference, 
A = m„ — m. p = 1.2935 MeV. It turns out that neutron- 
rich ejecta are only obtained when e Pe — e Ue > 4 A [TTI.[l2l| 
in the presence of similar electron neutrino and antineu- 
trino luminosities. 

The treatment of neutrino transport and neutrino 
matter interactions improved over the years signifi- 
cantly. Simulations of core-collapse supernovae, that in- 
clude three- flavor Boltzmann neutrino transport, showed 
rather small neutrino energy differences between D e and 
v e and hence generally proton-rich conditions were ob- 
tained during the early evolution after the onset of ex- 
plosion [22|, [23[ . More recently it has been possible to 
simulate core-collapse supernovae consistently through 
all phases up to 20 seconds after the onset of explo- 
sion [24], [25| , using general relativistic radiation hydro- 
dynamics that employs three-flavor Boltzmann neutrino 
transport and a sophisticated equation of state. It has 
been shown that the energy difference between v e and 
v e decreases continuously after the explosion has been 
launched and hence generally proton-rich conditions re- 
main for more than 10 seconds. This rules out initial 
expectations that the ejecta can become neutron-rich at 
late times. It leaves the fp-process [26|-[28j as the only 
nucleosynthesis process that can occur in neutrino-driven 
winds from non-rotating and not magnetized PNSs. 

In this article, we analyze the neutrino spectra ob- 
tained in the current simulations focussing in understand- 
ing the reasons behind the reduction of the energy dif- 
ference between v e and v e during the deleptonization of 
the PNS. We perform a detailed analysis of the main 
reactions determining the opacities or inverse mean free 
paths for the different neutrino flavors. Special attention 
is devoted to neutrino decoupling from matter, in par- 
ticular the determination of the position of the neutri- 
nospheres which evolves during the PNS deleptonization 
to continuously higher densities. We also present an al- 
ternative analysis, based on the formalism introduced in 
ref. that determines the contributions to the neu- 
trino luminosity from different processes and regions of 
the PNS. Both approaches show that at late times dur- 
ing the deleptonization, the spectra of all neutrino species 
forms at very similar radii and is determined by neutral- 
current processes that are insensitive to the neutrino fla- 
vor. 

The paper is organizes as follows. In section II we 
will introduce the procedure used to calculate the in- 



verse mean free paths for the relevant weak interaction 
processes considered, based on the neutrino distribution 
functions obtained in three-flavor Boltzmann neutrino 
transport. We will furthermore discuss neutrino decou- 
pling and introduce the concept of a scattering atmo- 
sphere following the work of ref. 30]. In section III, we 
will apply this analysis and discuss the evolutionary be- 
havior of the neutrino spectra obtained in core collapse 
supernova simulations of the low-mass 8.8 M Q O-Ne-Mg- 
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24j. In section IV we extend the analysis from sec- 
tion III to the 18 M Q iron-core progenitor from ref [24| 
and find qualitatively similar results. In section V is show 
the evolution of the neutrino spectra at different times 
during the deleptonization. While we explored opaci- 
ties in section III and IV, in section VI we analyze the 
neutrino emission, focusing on the location at which dif- 
ferent contributions to the luminosity at infinity appear. 
We close the manuscript with a summary in section VII. 



II. NEUTRINO SPECTRA 

We analyze data obtained in the core collapse super- 
nova simulations of massive stars published in ref. [24j. 
Stars in the mass range of 8.8-18 M Q were evolved con- 
sistently through core collapse, bounce and explosion in- 
cluding the neutrino-driven wind phase for more than 
20 seconds after core bounce. It includes most part of 
the PNS deleptonization, i.e. cooling via the emission 
of neutrinos of all flavors. For these models, the bary- 
onic equation of state from Shen et al. was used [3l|. 
It is based on relativistic mean field theory and the 
Thomas-Fermi approximation for heavy nuclei. Contri- 
butions from (e - , e + ) and photons where added following 
ref. H3- 

The simulations were based on general relativistic ra- 
diation hydrodynamics and three-flavor Boltzmann neu- 
trino transport in spherical symmetry [29l 1331-35] . The 
model is based on the following line element, 



ds 2 



-a 2 dt 2 



da 1 + dfl 



(1) 



which describes non-stationary and spherically symmet- 
ric spacetime, in coordinates t (system time) and a (en- 
closed baryon mass) , where dfl = r 2 (d9 + sin 2 9d(j)) de- 



scribes a 2-sphere of radius r(t,a). a(t,a) and T(t,a) 
are the metric functions and r(t,a)' — dr/da. With the 
following stress-energy tensor, 



■J), 



T tt = p(l -i 

T aa =p + pK, 
T ee = T^ =P+^p{J - K), 



(2a) 
(2b) 
(2c) 

(2d) 



with rest-mass density p, internal energy density e 
and matter pressure p, the equations for energy- and 
momentum-conservation are obtained via the covariant 
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derivative of the stress-energy tensor as follows, Vfc T kl = 
0. The quantities J, H and K are the zero, first and sec- 
ond angular moments of the neutrino distribution func- 
tions, f u (t,a, p,E), for more details see (29[. The mater 
velocity is given by the following relation, u — dr/adt. 

Neutrino transport determines the evolution of these 
neutrino distribution functions. In spherical symme- 
try, their momentum-dependency is given by the co- 
sine of the neutrino phase-space propagation angle (p) 
and the neutrino energy (E), for each neutrino flavor 
v 6 {v e , v & , Vfi/ri V/j,/t}- Their evolution is determined 
by solving the Boltzmann transport equation, assuming 
ultra-relativistic massless fermions, in a co-moving refer- 
ence frame. Defining the specific distribution functions, 
F v = f v /p, the Boltzmann transport equation can be 
expressed as follows [29|, |36| 



dF(p,E) 
adt 



H^-(Anr 2 apF(p,E)) 
aoa x ' 



1 1 da \ d r i o\ „, „,i 



(3a) 



(3b) 



^ + (30 



p a dr E 2 dE v " 



(3d) 



d In p 3u 



iA (£3jP(/i,jB)) (3e) 



adt 



(3f) 



collision 



Note that the neutrino distribution functions depend 
also on the system time t and radial coordinate (or en- 
closed baryon mass a), which are not shown here for sim- 
plicity. Expressions (j3aj) - (j3"e|) are the transport terms of 
the Boltzmann equation, while expression (j3f| is the col- 
lisional term. The latter is given by weak processes that 
change the neutrino distributions, as follows 



dF„(fi,E) 
adt 



collision 



\{E) 



E v (p,E) 



(4a) 



E 2 



c(hc) 3 
1 
P 



J dp' J d<t> R^ vN/A (E,E,cos6)F v (p',E) 



E 2 F„(p,E) 
c(hc) 3 



J dp' J d(jy R% tVN/A {E, E, cos 



- -FAp.E) ) / E' 2 dE' I dp' I d<f>R i » e ±(E,E',cos9)F u (p',E') 



F v {p,E)-^ J E' 2 dE' J dp' J d<f>R™l(E,E' , cos 9) (^-F v (n',E 

^3 J E' 2 dE' j 'dp' J # K,(E,E', cos9) Q-Fp (//,£') 
F„(ii,E)—^ J E' 2 dE' J dp' J dcj> R%(E,E',cos6)F p (p',E'). 



0) (4b) 



(4c) 



(4d) 



In the Boltzmann transport representation, weak pro- 
cesses are characterized by reaction rates. The relevant 
processes are listed in Table [I] For charge-current reac- 
tions these are the emissivity, j(E) and the inverse mean 
free path, 1/\{E). For neutral current processes we em- 
ploy scattering kernels, R{E, E' , cos 9). They define the 
probability of a neutrino of energy E and momentum 
space coordinates (p,, ip) being scattered to an energy 
E' and coordinates (p',ip'). Apart of the energies the 
scattering probability is a function of the scattering an- 
gle 13, 

COS0 = pp + — M 2 )(l — P' 2 ) COS0, (j) = if — if'. 

The elastic processes considered are neutrino scattering 
on nucleons (N) and nuclei (A), R\s.uN/a{E, E' , cos#). 



The inelastic processes are neutrino scattering on elec- 
trons/positrons, R™± ut (E, E' , cos 9), for in-scattering 

and out-scattering, and pair processes, R^J'{E, E' , cos 9), 
for the absorption of a neutrino pair (a) and for the pro- 
duction of a neutrino pair (p). The reaction rates for all 
processes will be further discussed below in section II A. 



A. Reaction rates and opacities 

Let us first focus on charge-current reactions and dis- 
cuss neutral current reactions further below. Assuming 
a composition given by free protons and neutrons, the 
change in the net electron fraction, Y e = Y e - — Y e + is 
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TABLE I. Neutrino reactions considered, including reactions rates and references. 





Reaction a 


Mean free path/Reaction kernel 


References 


1 


u e + n - 


■> p + e 




[37] 


2 


Ue+P ~ 


■> n + e + 




[3JJ 


3 


v e + {A,Z- 1) 






[3JJ 


4 


v + N - 


■> i/ 4- iV 


Rh, vN (E,E,cose) 


[3JJ 


5 


v + (A, Z) - 


+ 1/ + (4, Z) 


R^ A (E,E,co S 9) 


[3JJ 


6 


v + e ± - 


+ v' + e ± 


R™ ± (E,E', cos 9) 


[34, 35, 37] 


7 


v + v — > 


e~ + e+ 


Rtp {E,E', cos 9) 


[37, 38] 


8 


v + v + N + N^N + N 


Rt PNN (E,E', cos 6) 


[39] 



a Note: u = {ue,v e , f M / T , Vp/ T } and N = {n,p} 



given by 

Y e = - (r e - p + r PeP ) Y p + (r e+n + r„ c „) Y n . (6) 

where Y p and Y n are the proton and neutron abundances 
and the r's are rates for electron capture (e~p), antineu- 
trino absorption (v e p), positron capture (e + n) and neu- 
trino absorption (y e n). 

Boltzmann neutrino transport uses energy dependent 
neutrino emissivities, j{E), and opacities, x(E) (see 
equation (I4al0 . The latter is equal to the inverse mean 
free path, 1 / A(-E) . They depend on the distribution func- 
tions for and nucleons that are assumed to be in ther- 
mal equilibrium. The rates which appear in equation ([5]) 
are then expressed as follow 



2tt 



' e~p 



(hc) 3 n p 



dfi dE E 2 j Ve (E)[l- f Ve 0, E)] , (7a) 



?zl_ i dlldEE aj;>^m i (7b) 



(hc) 3 n 



(E) 



2tt 



r e +r 



(hc) 3 n 



- f d^dEE 2 j, c (E)[l~ fuA^E)}, (7c) 

n J 



2n 



(hc) 3 n n 



dfi dE E z c 



2 Ue(P> E ) 



A„ c (E) ' 



(7d) 



where n p and n n are the proton and neutron number 
densities. For the neutrino and antineutrino emissivities 
and mean-free paths we use the analytical expressions 
given in appendix C of ref. 3T~\. Note furthermore, the 
emissivity and opacity (or equivalent the mean-free path) 
are related by detailed balance as follows [33], 



ju e (E) = exp 



3v e ( E ) = ex P 



E -(fi e - hq) 
kT 



E+(ne- Mg) 
kT 



K B (E) ' 



(8b) 



where [1q — \i n — fi p . Note that we also use the reaction 
rates for charge-current processes with heavy nuclei from 
ref. [37]], which are based on a simplified description of 
the Gamow- Teller transition assuming an average nucleus 
with mean charge and mass. This has been improved 
recently, taking into account the distribution of heavy 
nuclei in nuclear statistical equilibrium [40L Elj . 

In a similar fashion as for the charge-current reactions, 
we can define mean free paths for the scattering pro- 
cesses: 



1 



\(E,n) c(h 



1 f°° r 1 

J dE'E' 2 j dv![l 
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R in (E,E',cost 



(9) 



They are given in terms of the reaction kernels for the 
different reactions, R(E,E',cos9) (see expressions ()4b[) - 
(UdJ). The scattering processes included in our simula- 
tions are scattering on nucleons and nuclei, electrons and 
positrons (reactions 4, 5 and 6 in table H]). Scattering on 
nucleons and nuclei is treated in the elastic approxima- 
tion also known as isoenergetic scattering (IS). In this 
case the scattering kernel reduces to R\s(E, E', cos 9) = 
R$ S (E,E, cos 9)5(E - E') and equation © gives: 



-E 2 



2 71 



djj! f dip 



Ms(E, n) c(hc)- 

xR? s (E,E,cos9)[l-f(n',E)}. (10) 

The expression we use for R® S (E, E, cos 9) is given in 
refs. [341 ■ For neutrino scattering on electrons and 
positrons (expression 6 in table HJ) we use the scattering 
kernels of references [HI IH, H3] ■ Note that in addition 
to having different kernels for v e and v e , we also include 
differences between v^/t an d v^ L / T . 

For pair processes (expressions 7 and 8 in table QJ the 
relation between the mean free path and the scattering 
kernel as follows 



1 



1 



\ v (E,fi) c(hc) 2 



dE'E' 2 / dn'f D (fj,',E') 



d<j)R%(E, E',cost 



(11) 



5 



and similarly for antineutrinos. Table U gives the refer- 
ences from where the expressions for the different scat- 
tering kernels have been obtained. Similar to equation 
the reaction kernels for in- and out-scattering as well as 
pair production (p) and pair absorption (a) are related 
via detailed balance [13] as follows 

iO (E, E> , cos 0) = exp | ~ {E kT E ' ] } Kf± (E, E', cos 0) 
R a vB (E, E', cos 0) = exp j E -^- J RF vB (E,E', cos 6) . 



B. Neutrino decoupling 

The neutrinosphere radii R V (E) at which neutrinos of 
energy E and flavor v decouple from matter, i.e. the 
surface of last scattering, is determined by the condition 
that the optical depth becomes 2/3. The optical depth, 
r, is defined as follows, 



Tv{E) 



n u (E) K(E)' 



(12) 



integrating the total inverse mean free path, 1/X V (E) = 
J2i l/Xi.v(E), from the stellar surface towards the center, 
where l/Xi }V refers to an individual weak process i. Fol- 
lowing refs [U we use the effective neutrino opacity 
for energy exchange or thermalization and define: 



1 



1 



A„,eff(#) 



Xv,abs(E) \A„ :a b s (-E0 A t / iSca tt(£') 



(13) 

The inverse absorption mean free path, 1 / Xi,,a.bs(E) , is 
obtained by summing contributions of processes in which 
neutrinos exchange energy with the medium (processes 
1-3 and 6-8 in table IJ) . The inverse scattering mean free 
path contains processes where momentum of the neutrino 
is changed but not their energy (processes 4 and 5 in 
table U). 

For those cases, scattering and pair processes, in which 
the opacity depends on both the energy and angle of the 
neutrino, we define the energy-dependent inverse mean 
free path as follows 



1 



X V (E) n v (E) 



dfj,f u (E,fi)- 



(14) 



X v (E,nY 

using the local neutrino distribution function for each 
neutrino flavor v and 



n v (E) = / dnf„{E,n). 



(15) 



In order to determine an ave rag e neutrinosphere ra- 
dius, the so-called energy sphere |30(, we define an energy 
averaged inverse mean free path using the local neutrino 
phase-space distribution functions: 



1 



1 2tt 

n v (he) 3 



J d i iE 2 dE f^E) 



1 

X v {E,nY 



(16) 
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FIG. 1. The region of neutrino decoupling is shown via the 
energy-dependent transport spheres (for z/ e ), where the color- 
coding is according to the neutrino energies from E m i n = 
3 MeV (black) at about 16.17 km up to _B max = 300 MeV 
(white) at 49.11 km. The energy- averaged z/ e -sphere is shown 
in the vertical blue solid line. The radial profiles of den- 
sity (solid black line) and luminosities (y^. solid red line, £> e : 
dashed red line, v^/t- dash-dotted line) are also shown. 



for each neutrino flavor v and 
2tt 



(he) 



diiE 2 dE f„(fJ,,E). 



(17) 



is the local number density of neutrinos. The energy 
averaged effective mean free path is defined as 



1 



Xu.cS 



1 



A,. 



abs 



l 



i/, abs 



l 



i>",scatt 



(18) 

Finally, the energy sphere radius i?gs is obtained from 
the condition: 



Tthcrm(-RES) = / dr / —— \ = \ 
I R ES \A„, c ff/ 3 



(19) 



For later discussions, it is convenient to define an aver- 
age transport neutrinosphere [3fl| Rtr from the condition: 



Ttr(Rtr) 



dr 



1 



(20) 
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where the total average transport mean free path is de- 
fined as follows 



Ai/,tr j \ ^iAabs / \ ^f,scatt / 

Note that if absorption processes dominate, 1/A^.tr > 
1/A„, c ff, while if scattering dominates 1/A„,tr >> l/A^cff, 
implying that the energy sphere is located as smaller 
radius (and hence higher density) than the transport 
sphere. As discussed by Raffelt [30], that means once 
neutrinos decouple form matter at the energysphere ra- 
dius, Res (eq.[19]), they still suffer several elastic scatter- 
ing events until they reach the transport sphere at larger 
radius R tr - The distance between both spheres defines a 
scattering atmosphere. 

Based on the above definitions, we will illustrate in the 
following the situation in the proto-neutron star (PNS) 
atmosphere. Fig.Q]shows radial profiles of matter density 
(solid black line) and neutrino luminosities (red lines), 
which correspond to a typical post-bounce time between 
1-2 seconds after the onset of explosion (depending on 
the progenitor model). The energy sphere for v e , see 
expression (|19[) is shown by a vertical blue line in the 
upper part of the plot. It is located near the PNS 
surface where the density drops over several orders of 
magnitude. Note that the density gradient at the PNS 
surface continues to steepen during the ongoing PNS 
deleptonization, which will be further discussed in sec- 
tion III. To illustrate the dependence of neutrino decou- 
pling on energy we also show the energy dependent neu- 
trinospheres, where color coding is according to the neu- 
trino energy between 3 MeV (black) to 300 MeV (white). 
Low energy neutrinos of 3-5 MeV decouple already at 
high densities (~ 10 14 g cm -3 ), while high-energy neu- 
trinos (100-300 MeV) decouple at rather low densities 
(~ 10 8 g cm~ 3 ). The luminosities can vary until the 
final neutrino decoupling radius, outside which the neu- 
trino luminosities stay constant and hence the spectra 
are frozen. 

Note that the luminosity hierarchy shown in Fig.Q]and 
its detailed evolution over time depends on the neutrino 
opacities used (see for example |25|), however, as the 
proto-neutron star deleptonizes one observes that both 
the luminosities and average energies of all neutrino fla- 
vors converge to almost indistinguishable values. It is 
one of the objectives of this manuscript to address the 
physical origin of this behavior. 

III. O-NE-MG-CORE SUPERNOVA EXPLOSION 

We begin considering the 8.8 M O-Ne-Mg-core col- 
lapse simulations. The model was explored consis- 
tently though core collapse, bounce, explosion and proto- 
neutron star (PNS) deleptonization up to 7 seconds after 
core bounce [24j. The early onset of explosion for this 
model, at already 35 ms post bounce, is related to the 
special structure of the progenitor that lacks an extended 



high-density envelope. For details about the explosion 
mechanism, see refs. [24],[43[. We will start our analysis 
by summarizing the evolution of neutrino luminosities 
and mean energies. Further below we will examine im- 
portant aspects related to the proto-neutrons star (PNS) 
deleptonization. 



A. Proto-neutron star deleptonization 

Fig.[5]shows neutrino energy L v and number luminosi- 
ties L n ^ u as well as root-mean-square and mean energies, 
with respect to time after bounce (in log-scale). For a 
definition of these quantities, see ref. j44(. All plotted 
quantities are measured in the co-moving reference frame 
at a distance of 500 km from the center (see ref. 24]). 
The sharp luminosity and energy jumps at about 40 ms 
post bounce are related to Doppler shift effects, where 
due to the passage of the supernova shock, matter ve- 
locities suddenly change from infall to expansion. Dur- 
ing the early post-bounce phase until 40 ms, the elec- 
tron neutrino luminosities are dominated by mass accre- 
tion at the neutrinospheres. They reach values of about 
2 x 10 52 erg/s and stay relatively constant during the ac- 
cretion phase, after the deleptonization burst from the 
core bounce has been launched. Note that during the 
accretion phase, the electron neutrino energy and num- 
ber luminosities are larger than the electron antincutrino 
energy and number luminosities. However, during the 
early explosion phase the electron antineutrino energy 
luminosity becomes slightly larger than the electron neu- 
trino energy luminosity (see Fig. [5]). The following hi- 
erarchy holds Lv r > L Ve > L Ufi/r until about 350 ms 
post bounce, after which the electron flavor neutrino lu- 
minosities change ordering and the following hierarchy 
holds L v . > Lp c > L Pc . This hierarchy remains until 
the simulation is stopped at about 7 seconds post bounce. 
The luminosities of all flavors decrease by about one or- 
der of magnitude within the first second after the onset of 
explosion. For the number luminosities the following hi- 
erarchy holds L n ,v e > -^ra,i> e > L n ^ up to about 750 ms 
post bounce, after which L n>Ve > L n ^ v . > L n> p B . This 
continuous decrease of the neutrino luminosities reflects 
the ongoing proto-neutron star deleptonization. 

Of additional relevance for nucleosynthesis and flavor 
oscillation studies are the mean energies, shown in Fig. [5] 
(bottom) , in particular the difference between v e and v e . 
During the accretion phase before the onset of an explo- 
sion, (E v .) decrease slightly from 16 MeV to 15 MeV, 
(Ep e ) increases from 10 MeV to 12 MeV and {E„ a ) de- 
crease from 13 MeV to 8 MeV. After the onset of explo- 
sion, both (E Uc ) and (E p<! ) increase slightly from 7.5 MeV 
and 10 MeV to 8.5 MeV and 12 MeV at about 350 ms 
post bounce. Both mean energies increase and their dif- 
ference stays as large and even increases slightly. How- 
ever, the difference is not large enough to turn matter 
neutron rich. (E& , ) decrease continuously from about 
14 MeV at the onset of explosion to 8 MeV at 7 seconds 
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FIG. 2. Post-bounce evolution of neutrino energy and number 
luminosities as well as mean and root-mean-square (rms) en- 
ergies for the 8.8 Mq O-Ne-Mg-core supernova [2J|. We show 
explicitly the average energies for neutrinos and antineutrinos 
of all flavors. For the luminosities, we only show u e , i> e and 
v v,It because iw T cannot be distinguished from f M / T at the 
scale chosen. 



post bounce. After about 1 s post bounce, also (Ep e ) 
and (E Ue ) decrease continuously to 9 MeV and 8 MeV. 
The rms-energies are slightly larger than the mean ener- 
gies but follow the same behavior. The decreasing mean 
energies for all flavors indicates the ongoing deleptoniza- 
tion of the central PNS and hence cooling by neutrinos. 



Furthermore, the mean energies of all flavors become in- 
creasingly similar with respect to time during the PNS 
deleptonization (see Fig. [2]). 

The resulting evolution of the explosion and the neu- 
trino spectra is in good qualitative and quantitative 
agreement with the study of the Garching group of this 
low-mass progenitor model [25l |43| , applying a different 
equation of state and in addition a set of updated weak 
interactions processes. 

The evolution of radial profiles of selected quantities is 
illustrated in Fig [3] at several post-bounce times (1 sec- 
onds: red lines, 2 seconds: blue lines, 7 seconds: green 
lines). We focus on the radial domain near the neutri- 
nospheres (vertical dashed lines for v e and dash-dotted 
lines for v & in graph (b)), i.e. the region where neutrinos 
decouple from matter and where the far distance spec- 
tra are determined. The graphs (a), (c) and (d) show 
radial profiles of temperature, entropy per baryon and 
electron fraction, all of which decrease at the neutri- 
nospheres. This evolution is typical for the PNS delep- 
tonization and neutrino cooling including the slow proto- 
neutron star contraction. Note the rapidly rising electron 
fraction outside the neutrinospheres, which is related to 
the expansion of material in the neutrino-driven wind 
where Y e ~ 0.56 (see ref. [24| for a discussion). As the 
temperature reduces, the mean energy of neutrinos also 
decreases and the neutrinospheres move to higher den- 
sities and hence smaller radii, during the proto-neutron 
star deleptonization (see Fig.[3]graphs (a) and (b)), from 
R Ve = 22.19 km and Rn e = 21.51 km at 1 second post 
bounce to R Ve — 15.28 km and Rp e = 14.97 km at 7 sec- 
onds post bounce. 

In the following subsection, we will analyze the rea- 
son for the decreasing difference in the mean neutrino 
energies. 



B. Individual opacities 

Fig. 2] shows radial profiles of inverse mean free paths 
for the individual reactions considered, for v e (left panel) , 
v e (middle panel) and v^/t (right panel) at selected post- 
bounce times obtained during the PNS deleptonization. 

We will start analyzing the inverse mean free paths for 
(fi, r)-neutrinos (right panel in Fig. 2]). Note that they 
have no contributions from charge-current processes, 
they are only produced via the neutral-current pair- 
creation reactions (7) and (8) in Table HI The dominating 
inelastic contribution comes from iV-iV-Bremsstrahlung, 
only a tiny contribution comes from e _ -e + -annihilation, 
and scattering on electrons/positrons. All inelastic pro- 
cesses are smaller by several orders of magnitude than 
elastic scattering on neutrons (IS, z/ p / T n). Note fur- 
ther that elastic scattering on protons (IS, v^/rP) is also 
smaller than scattering on neutrons because protons are 
much less abundant than neutrons. 

The large scattering dominance implies that the 
transport opacity is greater than the effective opacity, 
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Radius [km] Radius [km] 

FIG. 3. Radial profiles for selected quantities, for the 8.8 Mq O-Ne-Mg-core [24|], at three different post-bounce times (1 second 
(red lines), 2 seconds (blue lines), 7 seconds (green lines)). The vertical lines in the graph (b) correspond to the position of the 
energy spheres for j/ e (dashed lines) and (dash-dotted lines). 



l/A,y M/Ti tr 3> I/A;/ , , e ff, see equations (fT8|) and (|2Tj).The 
corresponding ncutrinosphcres spheres are shown in 
Fig. [4] (energy sphere: vertical black dash-dotted line, 
transport sphere: vertical black dashed line). The sepa- 
ration of both defines the scattering atmosphere for v^j T , 
which is present already at early times after the onset of 
an explosion. 

The situation remains the same at later times for 
(/i, r)-neutrinos. However, the difference between 
1/Ais,^ /tP and l/A^^i increases with 1/A^ 



becoming as large as \j\ v 



while 1/A^ /T&(4/T con- 
tinuously decreases with times (see the right panels 
Figs.|4(a)}j4(c)l). 



The situation for v e is illustrated in the middle panels 
of Fig. |4j Here, in addition to the neutral-current in- 
elastic scattering and pair-creation processes, absorption 
on protons contributes to the energy exchange processes. 
However, their contribution is small and comparable to 
neutrino scattering on electrons/positrons at early times. 
The small value of l/A^p, is related to the low num- 
ber density of protons in the neutron-rich environment. 
At late times l/A PcP even decreases an order of magni- 



tude below 1/Ay , p , jvat. As the density at the energy 
sphere, Res, increases with time, the opacity for v e ab- 
sorption on protons decreases due to Pauli-blocking of 
final-state neutrons. This aspect is taken into account 
via nucleon-degeneracy factors in the charge-current re- 
actions rates [37] . 

As before for v^ir, the dominating process is here 
also v e scattering on neutrons. Hence similar to (/x/r)- 
neutrinos, l/A Pc ,t r > 1/A^eff already early after the 
onset of an explosion and a scattering atmosphere has 
developed. At early times, due to contributions from v e 
absorption on protons the energy sphere for v e is located 
at larger radius, i.e. at lower matter temperature, than 
the energy sphere for v^/ T . However, as the opacity for 
v e absorption on protons decreases with time the tem- 
peratures at the energy spheres of D e and ^ M / T become 
increasingly similar. This explains the convergence of 
mean energies and luminosities for D e and v^/ T at late 
times (see Fig. [5]). 

Note the constant opacity for v e absorption on protons 
outside the transport sphere at 7 seconds post bounce 
in Fig. 21 which is due to the constant number density 
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(a) 1 second after bounce 




(b) 2 seconds post bounce 




(c) 7 seconds after bounce 



FIG. 4. Inverse mean free paths for the individual reactions considered (IS, vN): isoenergetic neutrino- nucleon scattering, 
(v e n, v e p): charge-current reactions (comparing rates form ref. (37J (Bruenn) and ref. [3| (Reddy)), (vis): e~ e + - annihilation, 
(yvNN): ./V-jV-Bremsstrahlung) , for v £ (left panel), v £ (right panel) and v^/t (right panel), based on radial profiles and 
evolution of the 8.8 Mq O-Ne-Mg-core collapse supernova simulation [24J] . We show selected post-bounce times 1 s (top), 2 s 
(middle) and 7 s (bottom). Note the different density scale for the bottom panel. The vertical black dashed and dash-dotted 

U~~ rn ~t ^--^^^-^ mil ~r FTTH 
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of protons where the increasing abundance of protons 
is compensated by the decreasing density at the proto- 
neutron star surface. 



In contrast to v e and 



the situation is different for 



v e illustrated at the left panels in Fig. |4] The opacity 
for v e is initially dominated by charge current v e absorp- 
tion on neutrons, 1/A„ e „. It is greater than l/Ais,i/ e Ti by 
several orders of magnitude early after the onset of explo- 
sion at 1 second post bounce, see Fig. 4(a) As the total 



opacity is dominated by absorption processes, we are in a 
situation where 1/A„ ej t r ^ 1/A^ e , ff- Hence, a scattering 
atmosphere is not present for v e a few seconds after the 
onset of an explosion [3(j. However, around 7 seconds 
post bounce a scattering atmosphere starts to develop, 
see Fig. 4(c) Similar as for P e , the opacity for neutrino 
absorption on neutrons becomes smaller than the one for 
neutrino scattering on neutrons in the region around the 
neutrinospheres. It is due to Pauli blocking of final-state 
electrons. 

At the time when the simulations were stopped, the 
effective opacity of v e still contains substantial contribu- 
tions from charge-current processes. Hence the spectra 
of v e and v e , the latter contains negligible contributions 
from charge-current contributions, have not converged 
yet. As will be discussed in the next section, we expect 
full convergence of the neutrino spectra and luminosities 
for v e and v e at later times. 

Note that for previous analysis is based on charge- 
current rates from (37} computed in the so called elastic 
approximation, which assumes zero-momentum transfer. 
This assumption may not be accurate at high densities, 
where neutrons become degenerate and neutrinos have 
higher average energies. To check this issue, we have in- 
cluded the full energy and momentum-transfer dependent 
treatment of ref. [46] for the charge-current reactions. 
The corresponding opacity based on that treatment is 
shown in Fig. |4] (red dashed lines, labelled 'Reddy') for 
v e and v e , in comparison to the opacity form ref. [37j (red 
solid lines, labelled 'Bruenn'). At low densities, the dif- 
ferences between both opacities are minor. At high densi- 
ties the differences between both approaches increases at 
late times as matter becomes increasingly neutron rich 
and the temperature decreases, which in turn rises the 
neutron degeneracy. 

It indicates that Pauli-blocking of final states described 
above, relevant for the suppression of charge-current re- 
actions during the deleptonization, will even increase in 
the simulations when applying the improved treatment 
of ref. Hl|. 

Now, we discuss the evolution of selected quantities at 
the v e energy sphere illustrated in Fig.[5]during the delep- 
tonization. Panel (a) shows the evolution of the sphere 
radii for v e (solid line) and v e (dashed line) . The temper- 
ature evolution at the f e -sphere is shown in Fig. [5] (b) , 
which remains about constant at 5 MeV. Fig. [5] (c) and 
(d) show the evolution of electron fraction at the v e - 
sphere and also the v e and v e fractions at the corre- 
sponding spheres. The electron fraction drops rapidly 




FIG. 5. Post-bounce evolution of the energy spheres for v e 
and v e , electron and neutron chemical potentials as well as 
temperature and electron fraction at the v e -sphere. Further- 
more, the evolution of the v e and v e fractions is shown at the 
corresponding spheres. 



from 0.1 to 0.035 during the first second after the onset 
of explosion. At later times, the electron fraction evolves 
to saturation values of 0.03 (at Res for v e) and 0.025 
(at Res for v e ). The neutrino fractions are already very 
small inside the proto-neutron star at the onset of ex- 
plosion, with continuously decreasing values Y Ve ~ 10~ 2 
and Yp e ~ 4 x 10~ 3 at the corresponding spheres. Later 
during the deleptonization, the neutrino fractions contin- 
uously decrease to Y Ve ~ 3 x 10 -5 and Yp e ~ 2 x 10 -3 
at about 7 seconds post bounce. Furthermore, we also 
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show the evolution of the neutron chemical potential, fi n , 
at the i/ e -sphere in Fig. [5](b). It rises continuously and 
reaches positive values at about 3.5 seconds post bounce, 
after which the neutrons have Fermi energies greater than 
their restmass. ji n exceeds the temperature at about 
6 seconds post bounce, after which the neutrons become 
degenerate. Equally important is the evolution of the 
electron chemical potential, at the i/ e -sphere, which 
is shown in Fig. [5] (b) too. It rises from around 10 MeV 
to 60 MeV during deleptonization, while the tempera- 
ture decreases only slightly. The rising electron chemical 
potential together with the decrease of the mean neu- 
trino energy implies increasing importance of final-state 
electron blocking. 



IV. IRON-CORE PROGENITORS 

In this section, we will compare the results found above 
for the low-mass O-Ne-Mg-core collapse supernova sim- 
ulation with the more massive 18 M© iron-core progen- 
itor [45j]. Note that neutrino-driven explosions of stellar 
models other than the O-Ne-Mg-core cannot be obtained 
using the standard physics input in spherical symmetry. 
Hence, in order to trigger the explosion for the 18 M Q 
progenitor under investigation, neutrino absorption and 
emission rates were enhanced in the gain region. Here 
we focus on the cooling phase of the proto-neutron star 
after the onset of explosion, that is not affected by the 
enhanced rates used to trigger the explosion. We have 
checked this by doing reference simulations without en- 
hanced rates where differences obtained are less than 5%, 
compared to the results presented in this section. 

Note further that in contrast to the O-Ne-Mg-core 
model, which was evolved up to 7 seconds post bounce 
where the spectra of v e and v e have not fully converged 
yet, the 18 M Q model was evolved for more than 20 sec- 
onds post bounce. Here the spectra have further con- 
verged and are indistinguishable for practical applica- 
tions. 

Fig. [B] shows the evolution of the neutrino observables 
for the 18 Mq model, again measured in a co-moving 
frame of reference at a distance of 500 km. The sharp 
jumps in luminosities and mean energies at about 350 ms 
post bounce are due to the onset of the explosion, where 
mass accretion vanishes and matter velocities suddenly 
change from infall to expansion (compare with the O- 
Ne-Mg-core in Fig. [2] where the explosion is launched 
already at about 35 ms post bounce). The post bounce 
mass accretion phase lasts for several 100 ms for this iron- 
core progenitor. Consequently the central PNS is more 
compact than for the O-Ne-Mg-progenitor, with higher 
central density and temperature as well as lower proton- 
to-baryon ratio. Furthermore, the longer accretion phase 
and the more compact PNS lead to very similar v e and 
v e energy luminosities during the accretion phase (see 
Fig. O, with the same ordering L 9c > L„ e >• L v ,. 
The number luminosities follow the same ordering here 
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FIG. 6. Post-bounce evolution of neutrino energy and number 
luminosities as well as mean and rms-energies for the 18 Mq 
Fe-core progenitors from ref. [24| . 



as for the O-Ne-Mg core, L n ^ e > L n ^ e > L„ )i/(i . The 
opposite ordering holds for the mean energies, (E„ c ) < 
(£ Pc ) < (E V)i/T ) < (E Pit/T ). They slowly rise during 
the accretion phase from (E Va ) — 8 MeV to 10 MeV, 
(Eu e ) = 11 MeV to 12.5 MeV, (E^ i/t ) = 15.5 MeV to 
16 MeV at the onset of explosion. 

At about 350 ms post bounce, the luminosities decrease 
rapidly from 4 x 10 52 erg/s (i/ e , v e ) and 1.5 x 10 52 erg/s 
(^V/r) to 7 x 10 51 erg/s (for all flavors) at about 0.5 sec- 
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Radius [km] Radius [km] 

FIG. 7. Radial profiles of selected hydrodynamic variables at 1 seconds (red lines), at 5 seconds (blue lines) and 20 seconds 
(green lines) post bounce. The same configuration as Fig. [3] but for the 18 Mq progenitor model from ref. [241 ] . 



onds post bounce. They also change ordering L u . > 
Lp e ~ L Ve . The same holds for the number luminosi- 
ties, L n>Ve > L n ^ v . > L n< p e , At about 20 seconds post 
bounce, the energy and number luminosities converge to 
1 x 10 50 erg/s and 1 x 10 54 s _1 . The mean energies de- 
crease also rapidly shortly after the onset of explosion 
(see bottom panels in Fig. [5]). On timescales of several 
seconds, the mean energies of all flavors decrease contin- 
uously and become increasingly similar. This evolution is 
in qualitative agreement with evolution of the O-Ne-Mg- 
core collapse supernova, discussed in section III. How- 
ever, for the iron-core model under discussion here, the 
simulation was carried out for more than 20 seconds aftre 
bounce. This allows us to explore with additional details 
the spectral convergence of the different neutrino flavors 
at late times. 

The evolution of radial profiles of selected hydrody- 
namic quantities is illustrated in Fig. [7] at the three se- 
lected post-bounce times 1, 5 and 20 seconds for the 
18 Mq iron-core progenitor model. Illustrated are con- 
ditions near the neutrinospheres during the PNS delep- 
tonization. Similar to the 8.8 M Q model discussed in 
section III (see Fig. [3]) temperature, entropy per baryon 



and electron fraction shown in graphs (a), (c) and (d) 
decrease continuously with time, while density at the 
neutrinospheres in graph (b) increases. The v e and v e 
spheres move not only to higher density but also closer to- 
gether during the proto-neutron star deleptonization on 
timescales on the order of several seconds (see graph (b)). 
It indicates the same behavior as for the O-Ne-Mg-core. 

In the following we repeat the analysis from section III 
B and look at the individual opacities (inverse mean free 
paths) at selected times during the proto-neutron star 
deleptonization, illustrated in Fig[5J Identical as for the 
O-Ne-Mg core, the opacity for v^j r is dominated by scat- 
tering on neutrons at any time. Inelastic processes, which 
are more than one order of magnitude smaller than elas- 
tic scattering on neutrons, are dominated by scattering 
on electrons/positrons. Only at late times after about 
5 seconds post bounce, the opacity of Bremsstrahlung 
rises and becomes of equal importance. The positions 
of energy sphere (dominated by inelastic processes) and 
transport sphere (all processes contribute equally) are 
marked by the dash-dotted and dashed lines in Fig. [5] 
Their separation indicates the presence of a scattering at- 
mosphere, which is already present for f M / T and v e early 
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(a) 1 second after bounce 




(b) 5 seconds post bounce 




(c) 20 seconds after bounce 



FIG. 8. Individual inverse mean free paths for the weak processes considered, based on the 18 Mq iron-core model [24|], at 
selected post-bounce times (1 second: red, 5 seconds: blue, 20 seconds: green). Vertical dashed and dash-dotted lines show the 
positions of transport and energy spheres. The same configuration as Fig. [4] 
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after the onset of explosion. Similar to the O-Ne-Mg- 
core, the situation is different for v ei for which the domi- 
nating opacity is absorption on neutrons. Consequently, 
energy and transport spheres lay close to each other and 
a scattering atmosphere has not yet developed. However, 
during deleptonization, the opacity for v e absorption on 
neutrons reduces continuously due to final state electron 
blocking. At about 10 seconds after bounce, the opac- 
ity from v e scattering on neutrons becomes larger that 
the opacity for absorption on neutrons, and a scattering 
atmosphere starts to develop. At about 20 seconds post 
bounce, the opacity of all flavors is dominated by neutral- 
current processes that do not distinguish between differ- 
ent flavors. The energy and transport spheres converge 
for all flavors and their spectra become increasingly sim- 
ilar. The evolution of opacities (inverse mean free paths) 
is in qualitative and quantitative agreement with the evo- 
lution for the O-Ne-Mg core discussed in section III B. 
We also find that the improved treatment of [46| leads to 
an additional reduction of the charge-current opacities 
(see Fig.©. 

Note that the extension in density of the scattering at- 
mosphere, i.e. the density domain between energy and 
transport spheres, reduces during the proto-neutron star 
deleptonization. However, measured in radius, the scat- 
tering atmosphere increases in size. 



V. NEUTRINO EMISSION CHARACTERISTICS 

In the previous sections we analyzed the neutrino spec- 
tra based on the opacities, i.e. the propability for a neu- 
trino to be absorbed on its way out of the proto-neutron 
star. Here we focus on neutrino emission processes, i.e. 
the inverse of absorption processes, and their contribu- 
tion to the luminosities at infinity. We follow the formal- 
ism of ref. [29( where the authors analyzed the early post- 
bounce phase up to 500 ms. According to appendix B of 
ref. (2^, we define the quantity Ci(E,r) such that the 
luminosity at infinity is given by 



Lna — 



dE 



Ci(E, 



(22) 



where i runs over all processes that create a neutrino of 
energy E at radius r. The quantity Ci(E,r) contains 
information about the local production of a neutrino of 
energy E and its absorption during the transport r — > 
oo. 

The function Ci(E,r) for the different processes con- 
sidered is shown in the Figs. [9] (8.8 M Q model) and [10] 
(18 M Q model), at selected post-bounce times. The col- 
ors of the contour lines represent the different processes, 
including charge current (blue), scattering on (green) 
and pair emission (white). Charge-current processes are 
the inverse of reactions 1-3 of table Q] and pair processes 
are the inverse of reactions 7 and 8 of table HI Note that 
elastic processes are not considered as they do not change 



the energy. The color gradient reflects the neutrino ener- 
gies, which range from 3 to 300 McV. The vertical lines 
at the top of each panel mark positions of the energy 
spheres and the dashed black lines are the luminosity 
profiles (left-axis scale). Furthermore, the vertical black 
lines at the bottom of each panel mark positions of dif- 
ferent densities in a logarithmic scale. 

As was discussed in ref. (see Figs. (4)-(7)), charge- 
current processes dominate the emission of v e and v e dur- 
ing the entire accretion phase before an explosion was 
launched. Here we analyze the later evolution during the 
proto-neutron star deleptonization. 

For the O-Ne-Mg-core illustrated in Fig. |9l v^/ T are 
produced by pair proc esses and scattering on e . At 
early times shown Fig. 9(a)| v^j T are produced by pair 
processes and their energy is modified by scattering on 
e ± . Furthermore, pair processes occur at high densities 
and hence contribute to neutrinos of higher energies. The 
produced high-energy neutrinos are down-scattered on 
at lower energies at low densities. The situation for 
pair processes and scattering on is similar for P e , how- 
ever the dominating emission contribution to the lumi- 
nosity at infinity comes from positron captures on neu- 
trons. It shifts the main production site to lower den- 
sities (larger radii) than for Vfi/ T - For v e , the emission 
is completely dominated by electron captures on protons 
at even lower densities. The different production sites 
reflect the hierarchy of the average neutrino energies. 

At later times, illustrated in Fig. |9(b)| we see that 
charge-current processes for v e become rather small. 
Hence, v e and v^/ T are mainly produced via the same 
neutral-current processes and at similar densities with 
slightly different contributions from pair processes and 
scattering one ± . It explains the slight difference between 
the T> e and v ll / T spectra (see Fig. [2]). For 24, we observe 
two clearly separated production sites. Pair processes 
and scattering on occur at almost the same densities 



as for v„ and v, 



'ix/t- 



However, the second contribution oc- 



curs at lower density (larger radius) via electron captures 
on protons. Both contributions have similar magnitude 
with the latter one producing neutrinos of lower energy. 

For the more massive 18 M Q model, the situation is 
similar at 5 seconds post bounce (Fig. 10(a)). Note the 



additional small contribution from electron captures on 
heavy nuclei, at densities on the order of 10 11 g cm~ 3 , 
which are expected to be present at the neutron-star 
crust. At later times, illustrated in Fig. |10(b)| charge- 
current contributions for v e are negligible and the spec- 
tra for v e and kW T become increasingly similar with time. 
For v e , we observe that charge-current contributions de- 
crease continuously and will become also negligible if the 
simulation was carried out further. It explains the slight 
difference remaining between the average energies of v e 
and v e (see Fig.©. These finding are consistent with the 
previous discussion based on opacities. 

It is interesting to note that pair processes, mainly N- 
A-Bremsstrahlung, occur at high densities on the order 
of 10 14 g cm -3 , where the treatment of nuclear correla- 
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FIG. 9. Local contributions, d{E,r), to the luminosity at infinity (see text for its definition) at two different post-bounce 
times during the protoneutron star deleptonization 8.8 Mq. Note the different t/-scales for each neutrino flavor. The color 
of the contour lines mark different weak processes (blue: charged current, white: pair processes, green: ^e ± -scattering). The 
color gradient reflects the neutrino energy from 3 to 300 MeV (see color bar at the bottom). The dashed black lines show the 
luminosities and the vertical lines the position of the energy spheres (equation 13). The black vertical lines at the bottom of 
each panel mark different baryon densities, indicated by the labels from log 10 (p [g/cm 3 ]) = 9-14. 



tions may affect the neutrino opacities [49|. However, all 
neutrino flavors will be influenced equally without pro- 
ducing any changes in the relative energies between dif- 
ferent flavors. 



VI. SUMMARY AND CONCLUSIONS 

We have performed a detailed analysis of the differ- 
ent processes that determine the neutrino spectra of all 
flavors during the deleptonization phase after the onset 
of supernova explosion. We have explored the 8.8 M 
O-Ne-Mg-core and the 18 M Q iron-core progenitors from 
ref. 2J], in order to cover a broad range of stellar models. 

Using neutrino opacities (inverse mean free paths) for 



the different processes, we computed energy and trans- 
port spheres for each neutrino flavor. Our results con- 
firm the finding of ref. pll. |30| for z-^/r, i.e. the presence 
of a scattering atmosphere between the v^/r energy and 
transport spheres. In addition we also find a scatter- 
ing atmosphere for v e already early after the onset of 
an explosion, which has not been observed before. We 
find that at early times after the onset of an explosion, 
due to the contribution from charge-current reactions to 
the total opacity for z/ e , the energy spheres of v e and 
Vfi/r are located at different positions. However, due to 
Pauli blocking of final state neutrons the £? e -opacity from 
charge current-reactions decreases continuously and the 
location of the energy spheres for v e and v^/ T becomes 
increasingly similar. At late times, the dominating con- 
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(a) 5 seconds post bounce (b) 22 seconds post bounce 



FIG. 10. Local contributions, d(E,r), to the luminosity at infinity (see text for its definition) at two different post-bounce 
times during the protoneutron star deleptonization for the 18 Mq model. The same configuration as Fig. [9] 



tribution to the total opacity of P e and v^jr comes from 
neutral-current reactions which are independent of the 
neutrino flavor. It explains the increasing similarity of 
the v e and v^/t spectra. 

On the other hand, the total opacity for v e is domi- 
nated by charge-current neutrino absorption on neutrons 
at early times and hence no scattering atmospheres is 
present yet. As the charged current opacity for z/ e is 
larger than for v e , due to the much larger abundance 
of neutrons compared to protons, neutral-current pro- 
cesses become dominating only at later times than for 
v e . The reduction in the charged current opacity for v e 
is related to Pauli blocking of the final state electrons, 
which continuously increase their Fermi energy as the 
neutrinospheres move to higher densities. At the same 
time, the average neutrino energies decrease. At these 
late times, a scattering atmosphere also develops for v e , 
which has also not been observed before, and the spectra 
of u e and v e become increasingly similar. 



In addition to the opacities, we have also explored 
emission processes (inverse opacities) from which we find 
qualitatively a similar behavior. Furthermore, we find 
that charge-current and neutral-current processes con- 
tribute to the v e luminosity at clearly separated den- 
sity domains. Charge-current contributions decrease dur- 
ing the deleptonization and originate from low densities 
around 10 10 -10 13 g cm -3 . On the other hand, neutral- 
current contributions originate at high densities between 
10 13 to several times 10 14 g cm -3 (depending on the pro- 
genitor and the state of deleptonization). 

A clear extension of this investigation is the inclu- 
sion of improved neutrino interactions, in particular cor- 
rections from weak magnetism and nucleon recoil [47| . 
the emission of (/i, r)-neutrino pairs via the annihilation 
of trapped electron neutrino pairs and the reverse pro- 
cess [48| . These were included in the study of the O- 
Ne-Mg-core in ref. (25|, which was in qualitative agree- 
ment with our findings. In addition, the treatment of 
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nuclear correlations [49( and weak processes with light 
nuclei [13, [EU should also be improved. Nevertheless, we 
expect that none of these additions will change the find- 
ings discussed in this paper, i.e. on timescales on the 
order of tens of seconds charge-current reactions are sup- 
pressed due to Pauli blocking of final states and the neu- 
trino spectra are dominated by neutral-current processes. 
These processes do not distinguish between different neu- 
trino flavors and hence the neutrino spectra of all flavors 
become very similar at these late times. Moreover, we ex- 
pect that including the improved phase-space treatment 
of ref. [46] in long-term simulations of supernova explo- 
sion, will even enhance the suppression of charge-current 
processes. 

We expect that the neutrino-driven ejecta will always 
be proton rich. This excludes neutrino-driven winds from 
non-rotating and not magnetic proto-neutron stars as 
possible site for the production of heavy r-process ele- 
ments [Z > 56). Furthermore, the role of neutrino- flavor 
oscillations at late times for neutrino detection and nu- 
cleosynthesis becomes negligible in the presence of very 
similar neutrino spectra of all flavors. However at early 
times, when the spectra of the different flavors are still 



different, it may impact nucleosynthesis under proton- 
rich conditions [52], [53[ and neutrino detection. 
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